LOCAL STRONG SOLUTION TO THE COMPRESSIBLE 
MAGNETOHYDRODYNAMIC FLOW WITH LARGE DATA 



XIAOLI LI, NING SU, AND DEHUA WANG 

Abstract. The three-dimensional compressible magnetohydrodynamic (MHD) isen- 
tropic flow with zero magnetic diffusivity is studied. The vanishing magnetic difTusivity 
causes significant difficulties due to the loss of dissipation of the magnetic field. The exis- 
tence and uniqueness of local in time strong solution with large initial data is established. 
The strong solution has weaker regularity than the classical solution. A generalized Lax- 
Milgram theorem and a Schauder-TychonofT-type fixed point argument are applied with 
novel techniques and estimates for the strong solution. 



1. Introduction 

Magnetohydrodynamics (MHD) concerns the motion of conducting fluids, such as gases, 
in an electromagnetic field. If a conducting fluid moves in a magnetic field, electric fields 
are induced and an electric current flow is developed. The magnetic field exerts forces 
on these currents which considerably modify the hydrodynamic motion of the fluid. On 
the other hand, the development of electric currents yields a change in the magnetic field. 
There is a complex interaction between the magnetic field and fluid dynamic phenomena, 
and both hydrodynamic and electrodynamic effects have to be considered. The equations 
for compressible magnetohydrodynamics consist of the Euler equations of gas dynamics 
coupled with the Maxwell's equations of electromagnetic field. The applications of mag- 
netohydrodynamics cover a very wide range of physical areas from liquid metals to cosmic 
plasmas, for example, the intensely heated and ionized fluids in an electromagnetic field 
in astrophysics, geophysics, high-speed aerodynamics and plasma physics. 

The equations of three-dimensional compressible magnetohydrodynamic flow in the isen- 
tropic case have the following form ( [21 [221 [23] ) : 

Pt + V • (pu) = 0, (1.1a) 
(pu) t + V • (pu ® u) + VP = (V x H) x H + pAu + (A + /z)V(V • u), (1.1b) 
H 4 — V x (u x H) = 0, V-H = 0, (1.1c) 

where p = p(x, t) 6 IR + denotes the density, u = u(x, t) € M 3 the velocity field, H = 
H(x, t) 6 M 3 the magnetic field, and P(p) = Ap 1 the pressure with a constant A > and 
the adiabatic exponent 7 > 1. The viscosity coefficients p and A of the flow are constants 
satisfying p > and 2p + 3X > 0, which ensures that the operator — pAu — (A + p)V(V • u) 
is a strictly elliptic operator. The symbol (£> denotes the usual Kronecker tensor product. 



Date: January 5, 2011. 

2000 Mathematics Subject Classification. 35Q36, 35D05, 76W05. 

Key words and phrases. Magnetohydrodynamics (MHD), zero magnetic diffusivity, strong solutions, 
existence and uniqueness. 

1 



2 



XIAOLI LI, NING SU, AND DEHUA WANG 



Usually, we refer to (jl.la[) as the continuity equation (mass conservation equation), and 
(jl.lbp as the momentum conservation equation. It is well-known that the electromagnetic 
fields are governed by the Maxwell's equations. In magnetohydrodynamics, the displace- 
ment current can be neglected ( [22j [23]). As a consequence, (jl.lcp is called the induction 
equation. As for the constraint V • H = 0, it can be seen just as a restriction on the initial 
value of H since (V • H)t = 0. We remark that, the magnetic diffusivity in (jl.ip is zero, 
which arises in the physics regime with negligible electrical resistance, see [5]. 
We consider the Cauchy problem of (jl.ip with the initial condition: 

(p,u,H)(x,0) = (po,u o ,Ho)0r), x G R 3 , (1.2) 

and are interested in the existence of solutions to (|l.l|) - (|1.2p . When the magnetic diffusivity 
v ^ 0, (fTTTcjl is 

Hi - V x (u x H) = -V x (z/V x H), V • H = 0, 

and there have been many studies and rich results in the literature: see [3 HI El [H O Q31 
[T9 | ll8 j [2"3" l l33] and the references therein. When the magnetic diffusivity v = as in (jl.ip . 
the mathematical analysis becomes much more difficult due to the loss of dissipation of 
the magnetic field, and to our best knowledge there have been no results on existence 
of solutions (even in the incompressible case). The aim of this paper is to establish the 
local existence and uniqueness of strong solution to system (jl.ip with large initial data 
in the three-dimensional space R 3 . By a strong solution, we mean a triplet (p, u, H) 
with u(-,t) G W 2 > q (R 3 ) and (p(-,t),H(-,t)) G W 1,9 (R 3 ), 3 < q < 6 satisfying (H) almost 
everywhere with the initial condition (jl.2p . As for the global existence of classical solutions 
of the small perturbation near an equilibrium for compressible Navier-Stokes equations, we 
refer the interested reader to |28jl29j and the references cited therein. The global existence 
of strong solutions with small perturbations near an equilibrium for compressible Navier- 
Stokes equations was also discussed in [30j [32] . Also see the discussions and references in 
[30j [32] for other related results on strong solutions. 

Throughout this paper, the standard notations for Sobolev spaces W S ' P (R 3 ) (H S (R 3 ), 
when p = 2) will be used. For p G [l,+oo], we denote by L?(0, T;X) the set of Bochner 
measurable X— valued time dependent functions ip such that t h-> \\ip\\x belongs to L p (0, T), 
and the corresponding Lebesgue norm is denoted by || • ||_Lppn- Denote the Sobolev 
space W l >P(0,T;X) := {ip | ip G Z7(0,T;JT), p t G LP(0,T;X)}, and H l {<d,T-X) := 
Ty 1,2 (0, T; X). Precisely, we will establish the following result on existence and uniqueness 
in this paper: 

Theorem 1.1. Assume that 

po G W^iK 3 ) n H^R 3 ), u G (H 2 {R 3 )) 3 , H G (^(K 3 )) 3 n (tf 1 ^ 3 )) 3 

for some q G (3, 6], and 

a < po < fi, \\po\\ w i, qnH i + ||u ||#2 + HHolliyi^nH 1 < r o 

for some positive constants ct,/3, and r$. Then there exist positive constants T = T{tq), 
ai(T, ro, a), and Pi(T, ro, 0), such that the Cauchy problem (jl.ip - (jl.2p has a unique strong 
solution (p, u, H) on R 3 x (0, T) satisfying 

pGL°°(0,T; W 1 ' 9 (M 3 )n,H' 1 (M 3 )); pt G L°°(0,T;L q (R 3 )); 
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oti(T,r ,a) < p < /3i(T,r ,/3); 

uG (L 2 (0,T; W 2 > q (R 3 )nH 2 (R 3 ))) 3 ; u t G (L 2 (0,T; H^R 3 ))) 3 ; 

He (L°°(0,T; W 1 ' 9 (M 3 )n J ff 1 (]R 3 ))) 3 ; H t G (L°°(0, T; L q (R 3 ))) 3 . 

In addition to the difficulties due to the presence of the magnetic field and its interac- 
tion with the hydrodynamic motion in the MHD flow of large oscillation, another major 
difficulty in proving the existence is the lack of the dissipative estimates for the magnetic 
field and the gradient of the density. Since we are concerned with the strong solutions 
in W 2 ' q which have weaker regularity than the classical solution in H 3 , we need some 
new techniques and estimates to establish the existence. We first linearize (jl.ip . use the 
Lax-Milgram theorem to obtain the solution to the linearized system, then apply the 
Schauder-Tychonoff fixed point theorem to obtain the strong solution of (jl.ip with large 
data. 

This paper is organized as follows. Section 2, which is the main body of this paper, is 
devoted to proving the local existence of the system (jl.ip by Lax-Milgram theorem and a 
fixed-point argument. Section 3 will focus on the uniqueness of the solution obtained in 
Section 2. 



2. Local Existence 

We use the letter C to denote any constant that can be explicitly computed in terms 
of known quantities, and the exact value denoted by C may therefore change from line to 
line in a given computation. Similarly, e and 5 denote arbitrary positive constants, C e and 
Cs denote correspondingly positive constants depending on | and |. Set di = d/dxi, i = 
1, 2, 3 for x = (x\, X2, X3) G R 3 . If A = ((fly)) and B = ((6y)) are 3x3 matrices, then 

33 33 

A:B = j2J2 a v b v and H = (^ : ^) I = (EE4) 5 - 

1=1 j=l 1=1 j=l 

\u\ p = norm in the space L P (M 3 ) with p > 1; 
||u|| s = norm in the Sobolev space H S (R 3 ) of order s on L 2 (M 3 ); 
||u|| sp = norm in the Sobolev space W S,P (M. 3 ) of order s on L P (M. 3 ). 
For a given T > 0, denote 

Qt :=M 3 x [0,T]. 

Obviously, LP(0,T;LP(R 3 )) = LP(Q T ). 

In this section, we will prove the existence result in Theorem 11.11 For simplicity of 
notations, we set jx = ft = 1 and A = without loss of generality. The proof will proceed 
through four steps by combining a generalized Lax-Milgram theorem and a Schauder fixed 
point argument. To this end, we consider first an auxiliary problem. 

Set 



$ = {<j>(x,t) G R 3 I G (L 2 (0,T;ff 2 (M 3 ))) 3 , & G (L 2 (Q T )) 3 }, 
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with the natural norm And, for q G (3,6], define 

Hf = $n |#M) G M 3 | G (L 2 (0,T;VF 2 ' 9 (M 3 ))) 3 n (L°°(0, T; P 2 (R 3 ))) 3 } 
n {<KM) G R 3 | G (L 2 (0,T;P 2 (R 3 ))) 3 , «fc G (L°°(0, T; L 2 (M 3 ))) 3 , 
V0 t G (L 2 (Qt)) 9 , <A(0) =u }. 
Using the continuity equation (jl.laj) . the momentum equation (If .fbj) can be reduced to 
pu t + pu ■ Vu + VP = (V x H) x H + Au + V(V • u), 

where the notation u • Vu is understood to be (u • V)u. Meanwhile, under the constraint 
V • H = 0, the induction equation is equivalent to 

H 4 + u VH = (Vu - (V • u)l) H, 

where I stands for the 3x3 identity matrix. 

We need to find a triplet (p, u, H) that satisfies the following auxiliary problem 

Pi + V-(pu) = 0, (2.1a) 
pu t - Au - V(V • u) = -pu • Vu - VP + (V x H) x H, (2.1b) 
H< + u VH = (Vu - (V • u)l) H, V • H = 0, (2.1c) 
a.e. in Qt for any given T > 0, with the initial condition f)1.2|) such that u G 
p€L°° (0, T; W 1 '"^ 3 ) n P X (R 3 )) n P 1 (0, T; L q (R 3 ) n L 2 (R 3 )) , 
ai < p < fii (ai,/3i are suitable positive constants), 

and 

H G (L°°(0, T; ll^QR 3 ) n P^R 3 ))) 3 n (P^T; L 9 (R 3 ) n L 2 (R 3 ))) 3 . 

where u G and u G (P 2 (R 3 )) 3 , p G W rl ' 9 (R 3 ) n P 1 ^ 3 ) wi * h « < Po < 1, H G 
(Ty 1,9 (R 3 )) 3 n (P^R 3 )) 3 are given functions and q G (3,6]. 

2.1. Solvability of the density with a fixed velocity. Obviously, the existence of a 
unique solution p := p(u) of the continuity equation follows directly from the method of 
characteristics. Although this method requires that u G (C°(0,T;C 1 (R 3 ))) d and p Q G 
C^R 3 ), the estimates below in Lemma 12.11 hold under the above assumptions on u and 
po, i.e., u G \P, po G W 1,9 (R 3 ) n P X (R 3 ). So we use, for simplicity, a formal approach. 
The correct procedure is to consider a regularization of u and po , and then to pass to limit 
(similar to the argument for H below, and can also be found in Theorem 9.3 in [30] ) . 

Lemma 2.1. Under the same conditions as Theorem \l.l[ there is a unique strictly positive 
function 

p:=p(u) G P 1 (0,T;L 9 (R 3 )nL 2 (R 3 )) nr 3O (0,T;lU 1 ' <7 (R 3 )nP 1 (R 3 )) 
which satisfies (|2.1ap . Moreover, the density satisfies the following estimate: 

l|Vp|| L oo( jYL9(R3)nL 2 (iR3)) < (||/?o||yKM(R3)afP(ffi 3 ) + C||u||#vTj eX p fC||u[|vpVT] . (2.2) 
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(|2.1ap can be rewritten as 

o 
d 

Then, the explicit formula for p is 



dt 

X(t) = X 



j t p(t, X(t)) = -p(t, X(t))V • u(t, X(t)). 



J Q V-u(r,X(r)) dr) . 



p(i,x) = p (X(0))exp 
It follows that 

aexp^-J^ iV-u^X^Udr^ < p(t, x) < exp (j |V • u(r, ^(r))^^ . (2.3) 

Now applying the gradient operator V to ()2. la|) , and using 

V(Vp • u) = u • V(Vp) + Vu • Vp, 
where Vu • V p = (Vu) T Vp, we have 

Vp t + u • V(Vp) + Vu • Vp + pV(V • u) + V • uVp = 0. (2.4) 
Multiplying (|2.4p by |Vp| 9 ~ 2 Vp with 3 < q < 6 and integrating over R 3 , we get 

-^l V < + / (- a • ^1^1") + W'^P ■ ( Vu • V P) 
qdt J R3 \q 

+ p\Vp\ q - 2 Vp • V(V • u) + |Vp| 9 V -uj dx = 0. 
Bearing in mind, by divergence theorem, 

/ (u-V(|VpH + |Vp| 9 V-u) dx= [ V • (\Vp\ q u) dx = 0, 

JR3 J R 3 

we find, by Hoder's inequality, 

1 d -\v P \ q < ivh^CIpIooIvcv • q)|, + ivpyvuU) + |vh«v • a d* 

: <1 JS? 



< ivpir'dplooivcv • a)|, + ivpyvuu + |Vp|||v • a| 

< C|Vp|«|VuU + IVpl^lploolVCV • a) l9 



Since 

W 1 ' q {R 3 ) ^ L°°{R 3 ) as 3 < q < 6, (2.5) 

then we have 

^|Vp| 9 < (^IVplglVuU + |p|oo|V(V • U)|, 

< C7||u||2,9|V^| g + |/£>|oo|V(V-u)| ff , 

and, by Gronwall's inequality, 

|Vp| 3 < exp J ||u|| 2i9 dr)(|Vp | g + J |p|oo|V(V • u)| g dr^ . (2.6) 
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Finally, (|2.2p follows from (|2.6p and Hoder's inequality. The proof is complete. 

□ 

2.2. Solvability of the magnetic field with a fixed velocity. Due to the hyperbolic 
structure of ()2.1cp in terms of the magnetic field H := H(u) with the fixed velocity u, we 
can solve H through the following Lemma 12.21 

Let Aj(x, t), j = 1, n, be symmetric mx m matrices in M. n x (0, T), f(x, t) and Vq{x) 
two m-dimensional vector functions defined in W 1 x (0, T) and R n , respectively. For the 
Cauchy problem of the linear system in V G M m x (0, T): 

V t +Y,A j {x,t)d j V + B{x,t)V = f(x,t), 
k V(x,0) = V (x), 

we have 

Lemma 2.2. Assume that 

Aj G (C7(0,T;fl' s (lR n ))nC 1 (0,T;F s - 1 (M rt ))) mxm , j = l,...,n, 

B G (C(0,r;F s - 1 (]R n ))) mxm , / G (C(0, T; ff s (R"))) m , V G (# s (R n )) m , 
with s > S + 1 an integer. Then there exists a unique solution to (|2.7p . i.e, a function 

v g (c([o, t), # s (R n )) n c x ((o, t), # 5 - 1 (R n ))) m 

satisfying (|2.7p pointwise ( i.e. in the classical sense). 

Proof. This lemma is a direct consequence of Theorem 2.16 in [30] with Ao(a?,f) = I. □ 

Now, taking advantage of Lemma 12.21 we have 
Lemma 2.3. Under the same conditions as Theorem there is a unique function 

H :=H(u) G (ff 1 (0,T;L«(R 3 )nl/ 2 (R 3 ))) 3 n (L~(0,T;W rl ' 9 (R 3 )nfl' 1 (R 3 ))) 3 
which satisfies the equation ()2.1cp . Moreover, the magnetic field satisfies 

l|H||i / cx)(o,r;VK 1 ^(K 3 )nH 1 (M 3 )) < M|Ho||wi>9'(R3)ni7 1 (ffi 3 ) + C||u||*\/r^ exp(C||u||^ \/T). 

Proo/. Assume that u G (C 1 (0, T; C^°(IR 3 ))) 3 and H G (C£°(M 3 )) 3 . Rewriting (gjcj) 
in the component form (u = (u^ 1 ), u^ 2 ), u®), H = (HW, H^ 2 \ H^ 3 ))), employing the 
summation convention on the repeated indices, we have 

nf ] + q^Sj-hW = (a 7 -uW - (v • a)<%) h^, * = i, 2, 3. 

i.e, 

H t + (u^I) 0,-H + ((V • u)I - Vu) H = 0. 

Clearly, Aj(x,t) = u^ j \x,t)I, j = 1,2,3, B{x,t) = (V • u)I - Vu and f(x,t) = 
satisfy the assumptions in Lemma 12.21 then we get a unique solution 

oo 

Hep {(C 1 (0,T,i/ s - 1 (R 3 ))nC(0,T; J ff s (lR 3 ))) 3 }, 

s=3 
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which implies, by the Sobolev imbedding theorems, 

oo 

He fl (C 1 (0,T;C fc (K 3 ))) 3 = (C 1 (0, T; C°°(]R 3 ))) 3 . 
fc=i 

For u £ f , H £ (W^iiR 3 )) n (iT 1 (R 3 )) , by an argument of dense set, there are 
{u»}£=i C (C 1 (0,T;C7* , (K 3 ))) 3 J {H J£° =1 C (C °°(K 3 )) 3 respectively such that 

u n — )• u in vp, 

H 0n ^H in^CR 3 )) 3 ^ 1 ^ 3 )) 3 . 

Hence, 

u n ^uin (C([0,T]xB(0,a))) 3 , 
for any a > and .6(0, a) denotes the ball with radius a and centered at the origin. 
Similarly, there are {H n }^ =1 C (C 1 (0, T; C°°(1R 3 ))) 3 satisfying 

H nt + u n • VH n = ( Vu n - (V • u n )l) H n (2.8) 

with H n (0) = H 0n . Multiplying ([22]) by |H n | p - 2 H„ (p > 2), integrating over M 3 , by 
integration by parts, we obtain, 

-4h„|p = -- [ Un- V(\H n \v)dx + / |H n r 2 H n • (Vu n - (V • u n )l) Ji n dx 

Pdt P PJ R 3 J R 3 

= - f |H n | p V • u n dx + [ |H n | p ~ 2 H n • (Vu„ - (V • u n )l) H n dx 

V 7r3 J R 3 

< C\H n \ p p \Vu n U 

where H n • (Vu n — (V • u n )l) = H^(Vu„ — (V • u n )l). Then, by Gronwall's inequality 



and the imbedding (|2.5D . we get 



i.e., 



|H n |£ < exp (c J iViinloorfr) |H n (0)|£ < exp ((J J* ||u n || 2i9 dr^ |H 0n 
|H n | p < exp (c J ||u n || 2l9 dT^ |H 0n | p . 



Thus, by Holder's inequality, we have 

l|H n [U°?(LP(R3)) - exp(C [lOnll^VrJlHonlp. (2.9) 
In particular, if we choose p = q in (I2.9p . then 

l|H n ||L5?(L9(M3)) < exp(C ||u n ||^Vr)|H 0rl |g < OO, 

and, up to a subsequence, assuming that {H n }^_^ were chosen so that 

H n — >• H weak-* in (L°°(0, T; L 9 (M 3 ))) 3 . 
Moreover, letting p — > oo in (|2.9p . using the imbedding (|2.5p again, we obtain 

l|Hn||i«»(Q T ) < exp(C||u n ||$\/f ; )|Ho n |oo 

< Cexp(C||u n ||*VT)||H 0n || lig < oo. 
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Taking the gradient in both sides of (|2.8p . multiplying by | VH n | 9 2 VH n and integrating 
over IR 3 , we get, with the help of Holder's inequality and the imbedding (|2.5|) . 

>*■« 

= - / |VH n |"- 2 VH n : (V(u n • VH n ))dx 

Jr 3 

+ / |VH n | g - 2 VH n : (V((Vu n - (V • u n )I)H n ))dx 

JR 3 

3 

= - / |VH„f- 2 V (a,H n • di(u n • VH n ) - diH n ■ d t ((Vu n - (V • u n )I)H n ))dx 

3 

= - / I VH n |«- 2 V (9iHW (5 4 u n • VHW + u n • di(VH^)))dx 

3 ij=i 

3 

+ / |VH„r 2 V (^W^^-fV.ujy 

+ diB^-H^di{d k v^ - (V • u n )S jk ))dx 
<C [ \VU n \ q \Vu n \dx-- [ u n -V(\VU n \ q )dx + C [ \VU n \ q \Vu n \dx 

Jr 3 q J R 3 J R 3 

+ C [ IHjVHrJ^lVVu^cfa 

<C/ |VH n | 9 |Vu n |dx + C / iHnllVHnl'-^Wunlda; 

Jr 3 Jr 3 



(2.10) 



< CIVUnloolVHnl' + C|H n | 00 ||u n || 2 JVEy<f 1 

< c||u n || 2i(? |vH n |5 + ciitg^ivEq*- \ 

Using Gronwall's inequality, we conclude that 

|VH n | 9 <e X p(C^ ||u n || 2 , g dr) (jVH n (0)| 9 + C jf || 
and hence, 

|VH|, < liminf | VH n | g < exp((7 jf ||u|| 2 , 9 dr) (jVH | g + C ||u|| 2 , 9 drj . (2.11) 
Furthermore, 

ll H IUf?(wi.«(R3)) < exp(C||ii||*\/r) ( ||H ||i, 9 + C\\u\\q,VT) . 



Passing to the limit as n — > oo in (|2.8p . we show that (|2. lc|) holds at least in the 
sense of distributions. Therefore, H t € (L 2 (0, T; L 2 (M 3 ))) 3 , then H e (^(0, T; L q (R 3 ) n 
L 2 (M 3 ))) 3 . The proof is complete. □ 
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2.3. Local solvability of (|2,lb|) . Now we prove the existence of a solution of (|2.1b[) . 
First we consider the bilinear form E(u, 0) and the functional L(4>) defined by 

E(u,4>)= [ (pu t - Au-V(V-u), (j> t -k(A(j> + V(V ■ 0))) dt 
Jo 

-(u(O),A0(O)+V(V-0(O))), 

L{<j>) = - f (pu ■ Vu + VP - (V x H) x H, t - k(A(j) + V(V • </>))) dt 
Jo 

- (u o ,A0(O) + V(V-0(O))) 
with k = (WpWl^^Qt))" 1 f° r 4> £ ^- Here, and in what follows, (•,•) denotes the inner 
product in (L 2 (M 3 )) 3 . Obviously, L(4>) is linear continuous on <1> with respect to the norm 



$. Moreover, by the Cauchy-Schwarz inequality, we get 

P(0,0) = f (Iv^tli + fclA^ + VCV-^ll-fcOo^A^ + VCV.^))^ 
J o 

+ J (|V*(T)l! + |V0(O)H + |V • 0(^)11 + |V • 0(0)11) 



> 



J (Wtll + ^1 + v(v • 0)H - - |l A <A + v(v • 0)|^ ^ 



+ 5 (|V0(T)H + |V0(O)|| + |V • 0(T)|| + |V • 0(O)|) 
>C||0lll- 

Hence, by the Lax-Milgram theorem (see |15j). there exists a u € $ such that 

£(u,0) = L(0) (2.12) 

for every ^ 6 $. 

If we assume that is a solution of the problem 

f t -A;(A0 + V(V-0)) =0, 
1 0(0) = h(x) 

with /i(x) smooth enough, and replace in (|2.12p by 0, then we have 

(u(0) - u , Ah + V(V • h)) = 0, 

which implies u(0) = uo- Similarly, if is a solution of the problem 

f t -fc(A0 + V(V-0)) =g{x,t), 
\ 0(0) = 

with g smooth enough, replacing by in (|2.12p . then we get 

r T 

/ (pu t - Au - V(V • u) + pu • Vu + VP - (V x H) x H, g) dt = 0, 
Jo 

which implies that (u, p, H) satisfies (I2.ip a.e. in Qt- 

Next, we prove the higher regularity for u. To avoid tedious calculations and notation, 
we work directly with the derivatives with respect to t of u instead of its differential 



2 -J IVpu t ||dr + i|Vu|| + i|V-u|2 

<^|Vu(0)|i + i|V-u(0)|i + | sup /VlilVulldT 

z z 4 0<r<t JO 

+ f f PV • u t dxdr + [ (H • VH, u t ) dr + )- f [ \H\ 2 V ■ u t dxdr 
Jo it 3 Jo 2 j j R 3 

< i|Vu(0)H + i|V • u(0)|i + C sup |p|oo f |Vu|l|Au| 2 dr 

^ * 0<T<t JO 

+ Cs dr + J* |-^H|^|VH|1 dr + jf |H|^|H|^r 



(2.13) 
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quotients. Multiplying (|2. lbj) by u t , integrating over Q t and using the Cauchy-Schwarz 
inequality, we obtain 

J \y/pu t \l ^+^|Vu|| + i|V-u| 

= ^|Vu(0)|| + ||V • u(0)H + J ( - pu • Vu - VP + (V x H) x H, u t )dr 

< ^|Vu(0)H + i|V ■ u(0)H + f I (\\^ t \ 2 + -\^pu • Vu| 2 ) dxdr 
2 2 Jo JR3 6 4 

+ / / PV • u t dxdr + / (H • VH, u t ) dr + - f [ |H| 2 V • u t dxdr. 
Jo Jw. 3 Jo 2 J J R 3 

By Holder's inequality and the Gagliardo-Nirenberg-Sobolev inequalities 

|Vu| 3 < C|Vu|J|Au|J, |u| 6 < C|Vu| 2 , (2-14) 
it is easy to deduce from f)2 . 13|) that 



(2.15) 



+ 5 [ \y/pu t \l dr + 25 [ |V • u t \ 2 2 dr. 
Jo Jo 

Now we differentiate (|2. lb>[) with respect to t and get 
p t u t + pvL tt - Au ( - V(V • u t ) 

= -p t u ■ Vu - pu t Vu-pu- Vu t — VPt + H 4 • VH + H • VH t - V(H • H 4 ). 

Multiplying (|2.16p by u^, integrating over M 3 , bearing in mind the continuity equation 
(|2.1a|) and the induction equation (|2.1c|) . we obtain 
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1 d 1 f 

~-dVMl2 + o / Pt\"t\ 2 dx + |Vu t |l + |V • u t || 

< |p|oo|u|oo|Vu|2|Vu| 3 |ui| 6 + |jO|cjo|u|||V(Va)]2]ut|6 + \p\ oo | u| oo | u| 6 1 Vu| 3 1 Vu t | 2 
+ |\/p|oo|\/pUt|2|Vu|3|Ut| 6 + | V ^|oo|Vu t | 2 |u| 00 | v /pU*|2 + \\P' || c o \p t | 2 | Vu t 1 2 

+ (V x (u x H) • VH, u t ) + (H • VH ( , u t ) + |HU|H t | 2 |V • u t | 2 

= |p|oo|u|oo|Vu| 2 |Vu| 3 |u t | 6 + Iploolulll V(Vu)| 2 |ut|6 + |/0|oo|u|oo|u|6|VQ| 3 |VUi| 2 

+ IVPloo|\/M| 2 |Vu| 3 |u t | 6 + IV^IoolVUt^lQlooly^Uf^ + \\P'\\ C ° |ft| 2 |Vu t | 2 (2.17) 

loc 

- / H ■ (Vut)V x(uxH)i- / H t (Vu t ) H dx + IHyH^lV • u t | 2 

Jr 3 Jr 3 

< |/o|oo|u|oo|Vu| 2 |Vu| 3 |u t | 6 + Iploolulll V(Vu)| 2 |ut|6 + |/o|oo|u|oo|u| 6 |Vu| 3 |Vu t | 2 

+ |\/p|oo|\/P U t|2|Vu|3|ut| 6 + Iv^UIVUt^lulooly^OUt^ + \\P' \\c° o JPt k\ Vu ( | 2 

+ iHlooluUlVHlalVu^ + |H|^|Vu| 2 |Vu t | 2 + |H|^|V • u| 2 |Vu t | 2 
+ iHloolH^lVutla + |HUH t | 2 |V • u t | 2 . 

Integrating (12.17P with respect to t, taking advantage of the continuity equation and the 
Gagliardo-Nirenberg-Sobolev inequalities as (|2.14p . we find, since p G L°°(Qx), 



i r* 

^\Vp^t\ 2 2 + J (\Vu t \l + \V-u t \ 2 2 )dr 

< hVfWMtyll + C I IPlooduloolVulflAulllVutla + IVulilAulalVutlaJdr 
1 Jo 

+ C / IVploodv^Ui^lVuHIAulllVui^ + |u|oo|Vu t | 2 |^ut| 2 )dr 
Jo 

+ ! \\P'\\c? l^|2|Vu t | 2 dT+ / (iHUIuUlVHblVu^ + lH^lVublVu^ 

Jo loc Jo 

+ |H|^|V • u| 2 |Vu t | 2 + iHloolHtlalVutla + IHUH^IV • u t \ 2 )dr ( 2 .18) 
<^l\/p(0)ui(0)|2 + C 5 { sup A|u|L|Vu|||Au| 2 + |Vu|i|Au|l)dT 

+ sup IVpIL f lVP Q t|2|Vu| 2 |Au| 2 dr+ sup |^/p|L f |u|LlVP u tll dl ~ 

0<r<t Jo 0<r<t Jo 

+ f\\P'\\ 2 c? \Pt\Ur+ AlHlLlulLlVHll + lH^lVull + lH^lV-ull 

Jo loc Jo 

+ 2|H|^|H t |l) dT}+S6[ \Vu t \ 2 2 dr + 5 [ |V • u t \ 2 2 dr + 5 [ \Vu t \ 2 2 dr, 

Jo Jo Jo 

where 5 > is small enough. 

Summing (|2.15p and (|2.18p . and for suitable S (6 < yg), we can first obtain, by Gron- 
wall's inequality, 

\\VP u t\\L™(L 2 (R'i)) < C, 
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and secondly, 

\\^t\l + ^|Vu| 2 + ~|V • u|| + i J \^pu t \ 2 dr + ±J (|Vu t || + |V • U t |l)dT 

< ||v^(o)ii*(o)|i + ~|Vu(o)|| + ~|v • u(o)|| + c 5 jf |p|2 dr + c 8 J* \^n\l\vn\ 2 2 dT 

+ C sup Moo f |Vu|i|Au| 2 dr + C s { sup |p|*, / (|u|^|Vu|||Au| 2 + |Vu|||Au||)dT 

0<r<* Jo 0<r<t JO 

+ sup \y/p\lo(( IVpUi|2|Vu| 2 |Au| 2 (ir+ I |u|^|^/pui|2dr) + / ||P'||;l |pt||dr 

0<r<t Jo JO JO ' oc 

+ A|H|L|u|L|VH|^ + |H|i|Vu|| + |H|^|V ■ u|l + 2|H|L|H t |l + |H|L|H|2)dr} 

JO 



+ <5 / \Vu t \ldr, 



which implies 



^pu t G (L°°(0,T;L 2 (1 3 ))) 3 ; u t G (L 2 (0, T; ff 1 (IR 3 ))) : 
Vu G (^(O.TjL 2 '™ 3 ^ 9 



(2.19) 



On the other hand, rewrite (|2. lh>[) as 



- Au - V(V • u) = -pu t — pu ■ Vu — VP + H • VH — ^V(|H| 2 ) = v, (2.20) 

which is a system with a strongly elliptic left-hand side. 

First, by a classical result on elliptic systems, there exists a positive constant C such 
that 

||u|| 2j? < C|Au + V(V-u)| 9 . 

Then v G (L°°(0, T; L 2 (M 3 ))) 3 by (j2~T3|) and ([2~T9]) . thus u G (L°°(0, T; Pf 2 (M 3 ))) 3 . Sec- 
ond, since p is bounded from below and y/pu t G (L 2 (Qt)) 3 , we know that u t G (L 2 (Qj-)) 3 . 
Hence, by the Gagliardo-Nirenberg-Sobolev inequality, we get 

Hllz4(L9(R 3 )) < C||u t ||^ 2(QT) ||Vu t ||^ ( Q T) as q G (3,6], 

for some 6 G [0,1). This implies that, by (I2T9D . u t G (L 2 (0, T; L 9 (M 3 ))) 3 , then v G 

(L 2 (0,T;L q (R 3 ))) 3 . A similar consideration leads to u G (L 2 (0, T; W 2 ' q (R 3 ))) 3 . Hence, 
we can conclude that u G \P. 

2.4. Existence for (jl.ip . The above argument guarantees the existence and uniqueness 
of the solution to the system (12. If) which enable us to define the map u = G(u) given by 
the composition of 

/: u-)-p(u), g: u ->• H(u), d: (p(u), u, H(u)) ->• u. 

Obviously, the fixed point of G : \& — )• \l/ is the solution of the system (jl.ip . To find a fixed 
point of G, we will use the Schauder-Tychonoff fixed point theorem (Theorem 5.28, |31j). 
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Consider the set 

M = {(j){x,t) G M 3 | max(||0|| £ 2 (w a,, (R 8) nH 2( R s)j, |lv^^tlU??(i 2 (R 3 ))) 

H0llL°?(fra(R3))5 ll^tlLl, (iJi(R3))) <r} 

with 

r 2 = cr(|Au |i + |u |^|Vu |2 + [|po[|v^i,<,nHi + II H o H^i^niJ 1 )/« 
for some suitable sufficiently large constant a > 1, where q G (3,6] . Clearly, M is a 
compact set in (L 2 (Qt)) ■ As we are going to use a fixed point theorem, we need to show 
that G(M) C M and G is continuous in M with respect to the norm in (L 2 (Qt)) 3 . 

We first prove that G{M) C M for some T = T. Indeed, assuming u 6 M, from (I2.3h . 
(12. 6p . and fj2 . 1 1 f) . by Hoder's inequality, we know that for 3 < q < 6, 

aexp(— Cry/t) < p(x,t) < exp(Cr\/i), 
< |Vp|, < (|Vpo|, + exp (CrVi) r^j exp (CV^) , (2.21) 
k |VH|g < (| VH |, + CrVt) exp (Cry 7 *) . 

Hence, from ([2TT8]) and (p£2l) , it follows that 

1 1 /"* 

-\VpU t \ 2 2 + - J \Vu t \ 2 2 dT 

< ^l\/K0)u t (0)|| + C 5 t(r 2 + V^t|li»(L»(R3)) + C S t(r 4 + r 6 ) + Ctr 2 . 
Since if we multiply (I2.1bj) by u^, integrate over R 3 and let t = 0, then 



(2.22) 



\VmMo)\i 

= [ (Au + V(V • u ) - p u • Vu - VP(po) + H • VH - ^V(|Ho| 2 )) • u t (0)da 

JR3 ^ 



< -^=(|Au | 2 + |V(V • U )| 2 + |/Oo|oo|u |oo|Vu |2 



+ \\P'\\ c o |Vp |2 + 2|Ho|oo|VH | 2 )|Vp(0)u t (0)| 2 , 

loc ' 



i.e., 



| v^(0)u t (0)|2 <-^(|Au | 2 + |V(V • uo)|a + |poUuoUVuo| 2 

\JOL 

+ ||P'|U |Vpo|2 + 2|Ho|oo|VH | 2 ). 

loc ' 

Combining (I2.22P and (I2.23p . taking 5 and T suitably small, we derive that 



(2.23) 



ll^ U 'llL S (L 2 (R3))+ll Vu *lli 2 (Q T )^r 2 ' 



T 



2 2 ^ 2 

\\VP U t\\L^(L 2 (R3)) + ll U *ll£2.(Hi(K3)) < g r • {2.24) 
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Then, we estimate the norm ||u|| L oc^2( R 3)) and the norm ||u|| L 2_( H /2, 9 ( K 3)). Indeed, taking 
advantage of (|2.20p . on the one hand, we get 

||AU + V(V-U)|| L =0 (L 2 (R 3)) 

+ C|I v / ? IIl^(l 2 (r3)) + C|l v H|| L ^( i2 ( R 3)), 
which leads to 

II u IIl^(#2( R 3)) < r 2 

from a classical result for elliptic systems. On the other hand, we have, also by the classical 
result on elliptic systems, 

F \\u\\l q dt<C C Woe | Vug + |Vp|3 + |u t g + |H&|VH|J)dt 
Jo Jo 

< Cr A T + r 2 T 

< r 2 

for some suitable small T. Hence, we have shown that G{M) C M. 

Next, we prove the continuity of G in M. We observe that if {u n }^ =1 C M, then there 
exists a subsequence (still denoted by {u n }^Li) such that 

u„, — > u strongly in M as n — > oo. 

Let p n and p be the solutions of 

p nt + V • {p n VLn) = 0, p n (0) = p , 

and 

Pt + v • (pa) = o, p(o) = po, 

respectively. Denote p n = p n — p, then p n satisfies 

p nt + u n -Vp n + (u n -u)-Vp + p n V-u n + pV-(u n -u) = 0, p n (0) = 0. (2.25) 

Multiplying (|2.25p by p n , integrating over Qt, and applying GronwalFs inequality, it is 
easy to derive that 

|p n || < exp(CrT) / (|(u - u n ) • Vp|^ + |pV • (u n - u)||) dt, 
Jo 

which implies that p n — > p strongly in L°° (0, T; L 2 (]R 3 )) . 

Similarly, we can show that H n -> H strongly in (L°°(0, T; L 2 (M 3 ))) . 

Now let u ra and u be the solutions of (I2.1bp corresponding to u n and u with 

Un(0) = u(0) = u . 

Then we have, denoting U n = u n — u and U n = u n — u, 

p„U nt - AU n - VV • U n 

= -p n u 4 - p n U n • Vu n - p n u • Vu„ - pu • VU n (2.26) 
+ H n • VH n - H VH - VP(p n ) + VP(p). 
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Multiplying (|2.26p by U nt , integrating over Q^, and take advantage of the convergence of 
p n and H n , we can prove as a routine matter that 

VU„ -> strongly in (L 2 (Q T )) 9 

and 

yfp^Unt -> strongly in (L 2 (Q T )) . 
Due to the convergence of p n , we deduce that 

U nt -> strongly in (L 2 (Q T )) 3 . 
Hence, by using the identity U ra (i) = J Q * XJ nt dT (U n (0) = 0), we get 

U n ->■ strongly in (L 2 (Q 7 )) 3 . 

Thus, the map G is continuous in M. The existence of a local solution is completely 
proved. 



3. Uniqueness 

We proceed to prove the uniqueness of the solution by the same procedure as that used 
for the continuity of G. We have already proved that for 3 < q < 6, 

u t G (L 2 (0,T;L 2 (M 3 )nL 9 (]R 3 ))) 3 , 

Vp£ (L 2 (0,T;L 2 (M 3 )nL 9 (l 3 ))) 3 , 

VHG (L 2 (0,f;L 2 (M 3 )nI 5 (K 3 ))) 9 . 

Using the standard interpolation, we get 

Ut G (L 2 (0,T;L 3 (IR 3 ))) 3 , 

Vpe (L 2 (0,T;L 3 (R 3 ))) 3 , 

VH G (L 2 (0,T;L 3 (M 3 ))) 9 , 

where 

i _ e l-e 

3 ~ 2 + g 

Now, assume that ui, U2 satisfy (jl.ip for some T > and denote 

p := p(ui) - /o(u 2 ) = pi - p 2 , u := ui - u 2 , H := H(ui) - H(u 2 ) = Hi - H 2 . 
Then, we have 

p t + Vp • ui + Vp 2 • u + pV • ui + p 2 V • u = 0, (3.1) 
with p(0) = 0. Multiplying (|3.1|) by p and integrating over R 3 , we have 

ld_ 
2~dt ] 



Ml 



- I |p| 2 V-uidx+ / pVp 2 -udx+ [ |p| 2 V-uicix+ / pp 2 V ■ udx = 0. 

2 Jr3 7R3 J R 3 J R 3 
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Combining the Cauchy-Schwarz inequality, the Holder inequality and |u|6 < C|Vu 

get, 

j t \p\l < |V • uilooHl + 2|u| 6 |pVp 2 || + e|Vu|l + C £ |p 2 |LMl 

< |V • uxloolpll + e|Vu| 2 + C £ \ P V P2 \l + e|Vu| 2 + C £ |p 2 |^|p|^ 

5 

< |V • UxUlpll + e\Vu\ 2 2 + C £ \V P2 \l\p\l + e|Vu| 2 + C £ |p 2 |L|p|^ 
< m (e)\p\ 2 2 + 2e\Vu\l 

where 771(e) = |V • mU + C £ (|Vp 2 || + |p 2 |^), £ > 0. 
Similarly, we have 

H ( + ui VH + u. VH 2 = (Vui - (V • ui)l) H + (Vu - (V • u)l) H 2 , 

with H(0) = 0. Using the same technique as for p, we get 



2dt 



U\ 2 2 --f |H| 2 V-uidx+/ H-(u-VH 2 )dx 
/ H • (Vui - (V • m)l) H dx + [ H • (Vu - (V • u)l) H 2 dx, 



and 



^-|H| 2 < C\V ■ ui|oo|H|| + e|Vu| 2 + C e |H 2 |^|H|2 + e |Vu| 2 + C £ |H • VH 2 | 2 

< C|V ■ u.Unll + 2e|Vu|i + C £ |H 2 |^|H|i + C e |VH 2 |||H|| 
<7? 2 (e)|H|2 + 2e|Vu|i, 



where m (e) = C|V • uxU + C e (|H 2 |L + |VH 2 | 2 ), e > 0. 
For Uj, i = 1,2, 

f piu it - Aui - V ■ Ui = - Pi ui ■ - VP(pi) + ^ • VHi - ±V(|H;| 2 ), 
[ Uj(0) = u . 

It is easy to derive 

pill* - Au- V(V • u) 

= -pu 2i - piu • Vui - pu 2 • Vui - p 2 u 2 • Vu - VP(pi) + VP(p 2 ) 
+ Hi • VHi - H 2 • VH 2 - ^VdH!) 2 ) + ^V(|H 2 | 2 ), 

with u(0) = 0. 
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Using once more the same technique as for p, H, and bearing in mind the continuity 
equation, we deduce that 

^IVpi u I! + I Vu I! + |v-u|2 

= — - / |u| 2 V • (piui)dx - I (pu 2t ■ u + piu ■ Vui • u + pu 2 ■ Vui • u 

2 JR3 JR3 V 

+ p 2 u 2 ■ Vu • u + VP(pi) • u - VP(p 2 ) u - H VHi u - H 2 VH • u 
+ i(V(|H 1 | 2 )-V(|H 2 | 2 ))-u)dx 

= / Vu • u • (piUi)dx + / ( pu 2t ■ u — p\u ■ Vui • u — pu 2 • Vui • u — p 2 u 2 • Vu • u 

- VP(pi) ■ u + VP(p 2 ) u + H VHi u + H 2 VH u 

- VHi • H • u - VH • H 2 • u) dx 

< C7 e |pi|L|ui|L|u|2 + e\Vu\ 2 2 + e\Vu\l + C B |uat|I|p|l + IVPiUVuiUMl 

+ luaUlVuiloodpl 2 + |u| 2 ) + Ql^l^lual^lul 2 + e|Vu| 2 + C £ \\P'\\ 2 c0 \p\ 2 2 + e|Vu| 2 

loc 

+ e|Vu| 2 + C e |VHi|l|H|l + CedU^l + |H 2 |2j|H| 2 + s\V ■ u| 2 

< 5e|Vu| 2 + e|V • u| 2 + m(e)(\p\ 2 2 + |u|l + |H| 2 ), 

(3.6) 

where 

7fe(e) = C e (|/Ol|Ll u llL + \ u 2t\l + IVPlloolVuiloo + lualoolVuiloo + |p 2 |Ll u 2|L 

+ ll-f'llc' + I 1 1 + (Hil^ + (H^ 

loc 

Summing up (|3.2p , f|3.4p and (|3.6p , by taking £ < ^, we obtain 

|(IVPi«ll + \P\1 + |H||) + |Vu| 2 < 2(771(5) + 77,(5) + 773(5)) (|u| 2 + \ P \\ + |H| 2 ) 

<77(e,i)(IVPiu|l + Hi + |H||), 

where 

77(5,*) = 2(771(5) +772(5) + 773(5)) max{l, -exp(Crv^)}. 

a 

And, the integrability of 77(5, t) with respect to t on [0, T] comes from the regularity of 
Ui,u 2 and the estimates in Lemmas 12.11 and 12.21 for pi, Hj with i = 1,2. Hence, 

IV^ul 2 , + Ml + |H| 2 = for all t G (0, T) 

follows from GronwalFs inequality, and consequently, 

u = 0, p = 0, H = on Q T . 

The proof of uniqueness is complete. 
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